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Branching diffusion in inhomogeneous media 

L. Koralov* 



Abstract 



£\j \ We investigate the long-time evolution of branching diffusion processes (starting 

with a finite number of particles) in inhomogeneous media. The qualitative behavior 
of the processes depends on the intensity of the branching. In the super-critical 

O . regime, we describe the asymptotics of the number of particles in a given domain. 

In the sub-critical and critical regimes, we show that the limiting number of particles 
is finite and describe its distribution. 
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1 Introduction 



Consider a collection of particles in M. d that move diffusively and independently. Besides 
the diffusive motion, the particles can duplicate with the rate of duplication /3v(x), x e M. d , 
where x is the position of a given particle, v is a continuous non-negative compactly 
supported function and (3 > is a parameter controlling the duplication rate. Both 
copies start moving independently immediately after the duplication. 

There is a number of results describing the processes in the case when v in concentrated 
in one point. See, for example, [6], [8], [TJ, [9] for the study of superprocesses and PQ, 
!3] • • [IB] for the asymptotic properties of branching random walks. We now will give a 
detailed description of the behavior of the branching diffusion in M d in the case when v is 
a compactly supported function. Our main goal is to describe the distribution of particles 
when t is large. The asymptotics depends on whether j3 is above, at, or below the critical 
value /3 CT , which is the infimum of values of /3 for which the operator 

lfu(x) = -Au(x) + /3v(x)u(x) (1) 

has a positive eigenvalue. This is the operator in the right hand side of the equations on 
the particle density and higher order correlation functions, given below. 



*Dept of Mathematics, University of Maryland, College Park, MD 20742, koralov@math.umd.edu 



We will show that for j3 > (3 cr the number of particles in a given region U at time t 
has the asymptotics 

nt(U) ~ e A °^ / f(y)dy, 

where Ao(/3) is the largest eigenvalue of Zr, £ is a random variable that depends on the 
initial configuration of particles (assumed to be finite in number) and tp is a determin- 
istic function (limiting density profile). Intuitively, the presence of the random variable 
£ reflects the effect of branching at random times while the number of particles is small. 
After the number of particles becomes sufficiently large, it keeps growing nearly deter- 
ministically due to the fact that the bulk of particles is located near the support of v and 
due an 'averaging' effect in the branching mechanism. 

When (3 < (3 cr (and d > 3), the effect of transience of the diffusion outweighs the 
branching, and all the particles will eventually wander off to infinity (assuming that 
initially there were finitely many particles). It will be shown that for (3 < (3 cr the total 
number of particles tends, as t — > oo, to a finite random limit, whose distribution will be 
identified. 

The case when /3 = (3 CT is interesting when d > 3 (/3 CI = for d = 1,2). In this case the 
total number of particles will be shown to tend to a finite limit almost surely, although 
the expectation of the total number of particles tends to infinity. 

Some of the corresponding results for the processes on the lattice with branching at 
the origin were obtained in [I], [3], [2], [18] . The main difficulty here, compared to the 
latter series of papers, is that the explicit formulas for the resolvent of the generator that 
were helpful in analyzing the processes with branching at the origin are not available 
now. Besides allowing the treatment of the general potential v, the techniques developed 
in this paper will allow us to study some of the more intricate properties of the limiting 
distribution: the fluctuations of the local number of particles conditioned on the total 
number of particles in the super-critical case, the growth of the region containing the 
particles and the distribution of the number of particles in the regime of large deviations 
(near the edge of the region containing the particles). These properties will be the subject 
of a subsequent paper. 

The results on the large time asymptotics (Sections HJ[6]) will be obtained from the 
equations on the particle density and higher order correlation functions derived in Sec- 
tion [2j The analysis will be based on the spectral representation of solutions in the 
appropriate function spaces followed by the asymptotic analysis of integrals with inte- 
grands that depend on several parameters. Some of the techniques are related to those 
employed by us in the study of a polymer distribution in a mean field model ([5]). There, 
the asymptotics of a single equation (rather than a recursive system of equations) in the 
parameters t and /3 was examined. 

Finally, let us mention a number of recent papers on the parabolic Anderson model, 
where v is a stationary random field (see, for example, [T3], [TJ], [I], [E], [TT], [U])- When 
v is random, the behavior of the solution to ([I])-© essentially depends on nature of the 
tails of the distribution of v. It has been shown that in many cases the main contribution 



to the creation of the total number of particles is given by the isolated high peaks of 
the random potential. Moreover, when t is large, the bulk of the solution is located near 
one of those peaks with high probability. This adds to the importance of the study of 
branching diffusions in the case when v is localized. 

2 Equations on correlation functions 

Let B$ be a ball of radius 5 in M d . For t > and x, yi, 1/2, ... EM. d with all yi distinct, define 
the particle density pi(t, x, yi) and the higher order correlation functions p n (t, x, y±, ...,y n ) 
as the limits of probabilities of finding n distinct particles in Bs(yi),...,B$(y n ), respectively, 
divided by Vo\ n (B$), under the condition that there is a unique particle at t — located 
at x. We extend p n (t,x,yi, ...,y n ) by continuity to allow for y^ which are not necessarily 
distinct. For fixed y%, the density satisfies the equation 

d t pi(t,x,yi) = -Ap 1 (t,x,y 1 )+ /3v(x)p 1 (t,x,y 1 ), (2) 

Pi(0,x,y 1 ) = 6 yi (x). (3) 

Indeed, let s,t > 0. Then we can write 

Pi(s + i,x,yi) = (2-ks)~* I e^^p 1 {t,z,y 1 )dz + (3v{x)sp 1 {t,x,y 1 ) + a{s,t,x,y 1 ), (4) 



where the term with the integral on the right hand side is due to the effect of the diffusion 
on the interval [0, s], the second term is due to the probability of branching on [0, s], and 
a is the correction term. The correction term is present since (a) more than one instance 
of branching may occur before time s, and (b) even if a single branching occurs between 
the times and s, then the original particle will be located not at x but at a nearby 
point and the intensity of branching there is slightly different from (3v (x). It is clear that 
lim s ^ sup^^gjgd a(s, t, x,y)/s = 0. After subtracting p 1 (t,x,yi) from both sides of (j4j), 
dividing by s and taking the limit as s I 0, we obtain (T2J). 

The equations on p n , n > 1, are somewhat more complicated: 

d t p n (t, x, y u ..., y n ) = -Ap n (t, x, y h ..., y n )+/3v(x) (p n (t, x, y u ..., y n ) + H n (t, x, y ± , ..., y n )) , 

(5) 
p n (0,x,yi, ...,y n ) =0. (6) 

Here 

H n (t,x,y 1 ,...,y n )= ^ P\u\(t,x,U)p n _iui(t,x,Y\U), 

UCY,U^<D 

where Y = (yi, ..., y n ), U is a proper non-empty subsequence of Y, and \U\ is the number of 
elements in this subsequence. Equation (JSJ) is derived similarly to ([2]). The combinatorial 
term H n appears after taking into account the event that there is a single branching on 



the time interval [0, s], the descendants of the first particle are found at the points in U at 
time s + t, while the descendants of the second particle are found at the points of Y \ U, 
with the summation over all possible choices of U. 

3 Analytic Properties of the Resolvent 

Here we recall some basic facts about the operator L 13 : L 2 (W d ) — >■ L 2 (R d ) (see (CD)) and 
its resolvent R x = (LP — A) -1 . We will assume that v > is continuous, compactly 
supported and not identically equal to zero. It is well-known that the spectrum of LP 
consists of the absolutely continuous part (— oo, 0] and at most a finite number of non- 
negative eigenvalues: 

a(LP) = (-oo, 0] U {Xj}, 0<j<N, Xj = Xj(/3) > 0. 

We enumerate the eigenvalues in the decreasing order. Thus, if {A^} ^ 0, then A = 
maxAj. Thus the resolvent R x : L 2 (IR d ) —¥ L 2 (M d ) is a meromorphic operator valued 
function on C = C\(— oo, 0]. 

Lemma 3.1. There exists (3 cr > (which will be called the critical value of (3) such 
that supa(Z/) = for (3 < /3 cr and supa(L /3 ) = A (/3) > for /3 > /3 cr . For /3 > /3 cr 
the eigenvalue Xq(/3) is a strictly increasing and continuous function of (3. Moreover, 
lim^^ X(/3) = and lim^oo A(/3) = oo. 

The proof of this lemma is standard (see Lemma 4.1 of [5]). Denote the kernel of R x 
by R^(x, y). If (3 = 0, the kernel depends on the difference x — y and will intermittently 
use the notations Rf{(x,y) and -R°( x — v)- The kernel R\(x) can be expressed through 
the Hankel function Hi, . 

Rl(x) = c d k d - 2 (k\x\f-iH { p(iV2k\x\) ) k = ^X, Rek>0. (7) 

2 1 

We shall say that / G C exp (M d ) (or simply C exp ) if / is continuous and 

H/llc cxp (^) = sup(|/(x)|e N2 ) <oo. 

xeR d 

The space of bounded continuous functions on M d will be denoted by C(lR d ) or simply C. 
The following lemma will be proved in the Appendix. 

Lemma 3.2. The operator R x : C exp (IR d ) — > C(lR d ) is meromorphic in X E C. Lts poles 
are of the first order and are located at eigenvalues of the operator L^ . For each e > 
and some A = A(/3), the operator is uniformly bounded in X £ C , |argA| < it — e, \X\ > A. 
Lt is of order 0(1/ 1 A| ) as X — > oo, |argA| < n — e. 

Lf d > 3 and (3 G [0,/3 cr ), then R x has the following asymptotic behavior as X — > 0, 
X G C: 

^ = QjAi- 1 + P ( f(A) + 0(|A^), d>3, d-odd, (8) 



J R^ = Q2A^ 1 ln(l/A) + Pf(A) + 0(|A|5|lnA| + |A| d - 2 |lnA| 2 ), d > 4, d - even, (9) 

where P d are polynomials with coefficients that are bounded operators and Q d are bounded 
operators. 

The limit P^(Q) = limA->o,Aec R\ will be denoted by Rq. It is an operator acting from 

C cxp (R d ) to C(R d ). 

It will be shown in the Appendix that if /3 > /3 CT , then the eigenvalue Ao(/3) of the 
operator L^ is simple and the corresponding eigenfunction does not change sign (and so 
can be assumed to be positive). From this and Lemma [3.21 it follows that the residue of 
R x at A is the integral operator with the kernel ip p(x)ip p(y) , where ipp is the positive 
eigenfunction normalized by the condition ||^||L 2 (R d ) = 1- Note that tj)p decays exponen- 
tially at infinity. More precisely, it follows from (|7|) that if we write x as (9, \x\) in polar 
coordinates, then there is a continuous function fp such that 

ipp(x) ~ fp(9)\x\2~z exp(— y2Ao|x|) as |x| — >■ oo. (10) 

If /3 = (3 CT , then Ao = might not be an eigenvalue of the operator lr . As shown in 
Lemma T7.31 for d > 3 there is a unique (up to a multiplicative constant) positive function 
ipp (the ground state of L 13 ) which satisfies 

L% = -A^ + pvfo = 0, ippix) =0(\x\ 2 - d ), ?$2-(x) = 0(\x\ 1 - d ) as r = \x\ -> oo. 
2 or 

In fact, ipp is a genuine eigenvector (element of L 2 (M. d )) if and only if cf > 5. We will 
normalize i[)p by the condition that ||^f0/3||L 2 (K d ) — 1- 



4 The super-critical case 

Throughout this section we assume that (3 > (3 cr . First, let us introduce some notations. 
Note that the dependence of some of the quantities below on in not reflected in the 
notation in order to avoid overcrowded formulas. For a positive number x, we define the 
curve T(x) in the complex plane as follows: 

r(x) = {A : |ImA| = y/4x(x-ReX), ReA > 0} |J{A : |ImA| = 2x(l - ReA), ReA < 0}. 

Thus r(x) is a union of a piece of the parabola with the vertex in x that points in the 
direction of the negative real axis and two rays tangent to the parabola at the points 
it intersects the imaginary axis. The choice of the curve is somewhat arbitrary, yet the 
following properties of T(x) will be important: 

First, ReA < x for A G T(x). Second, since the rays form a positive angle with the 
negative real semi-axis, we have |argA| < n — e(x) for all A G T(x) for some e(x) > 0. 



Third, since the rays are tangent to the parabola, and the parabola is mapped into the 
line {A : Re A = yfx} by the mapping A — > a/A, the image of the curve T(x) under the 
same mapping lies in the half-plane {A : Re A > y/x}. 

The integration along the vertical lines in the complex plane and along contours T(x), 
below, is performed in the direction of the increasing complex part. 

We'll need estimates on the solutions of the following parabolic equation. Let 

d t p(t,x) = -Ap(t,x) + f3v(x)p(t,x), p(0,x) = g{x) e C exp . (11) 

We'll denote the Laplace transform of a function / by /, 

POO 

/(A) = (£/)(A) = / exp(-\t)f(t)dt. 



Jo 

Let r be the distance between Ao and the rest of the spectrum of the operator LP. In the 
arguments that follow we'll use the symbol A to denote constants that may differ from 
line to line. 

Lemma 4.1. For each e e (0, r), the solution of < f77]j has the form 

p(t, x) = exp(A t) (V>/9, g)ipp(x) + q(t, x), (12) 

where 

\\q(t,-)\\ c <A(e)ex V ((\ -6)t)\\g\\ Ccxp . 

Proof. After the Laplace transform, the equation becomes 

(-A + f3v)p-Xp = -g. 
Thus, the solution p can be represented as 

p(t,-) = ~[ e xt R{gd\. (13) 

A7[% </ReA=A +l 

The resolvent is meromorphic in the complex plane outside of the interval (— oo, Ao — r], 
with the only (simple) pole at Ao with the principal part of the Laurent expansion being 
the integral operator with the kernel il)p(x)ij}p(y)/(\Q — A). 

By Lemma |3.2[ the norm of the R x does not exceed A/|A| near infinity to the right 
of r(A — e). Therefore, the same integral as in (fl3"|) but along the segment parallel to 
the real axis connecting a point Ao + 1 + ib with the contour T(Ao — e) tends to zero 
when b — > cxo. Therefore, we can replace the contour of integration in f )13p by T(Ao — e). 
The residue gives the main term, while the integral over T(Ao — e) gives the remainder 
term. □ 



Lemma 4.2. Let K C M. d be a compact set. For each e G (0, r), the function pi(t,x,y) 
satisfies 

Pi(t, x, y) = e^q>(\ o t)ip0(x)ip0(y) + q(t, x, y), 

where 

sup \q(t,x,y)\ < A(e) exp((A - e)t - \y\ \/2(A - ej) 
xeK 

fort> 1/2. 

Proof. Let K' be a compact set that contains supp(t> ) U K in its interior. Consider first 
the case when y G K' . Apply ( 1T2|) with t replaced by t' = t — 1/2 and g = pi(l/2, -,y). 
In order to calculate the main term of the asymptotics, we note that ||g||c cx is bounded 
uniformly in y G K' and 

(^p,9) =exp(-A )^(l/)- 

The latter follows from 

/■1/2 Q 

= J % + ^)(exp(-A t)^), pi)dt = 

(exp(-A £)^,pi>|Eo + J ((exp(-Aot)M (-^ + L p ) Pl )dt = 

(exp(--A )^/3,Pi(l/2, •,?/)) - (^/9,pi(0,-,j/)) = 
1 

ex P(-2 A o)(^/3,^)-^/3(Z/)- 

Therefore, (fT2l) implies that 

Pi(t, x, y) = exp(A t)^,3(y)^/3(a;) + exp((A - e)t)q(t, x, y), 

where \\q(t, -,y)\\c < A(K') f° r a U V *= ^'- It remains to consider the case when y £ K' . 
Let u(t, x, y) = pi(t, x, y) —po(t, x, y), where p is the fundamental solution of the heat 
equation. Then u satisfies the non-homogeneous version of (fTTl) with the right hand side 
/ = —(3v(x)po(t,x,y) and g = 0. Note that / is a smooth function since y £ K' . Solving 
this equation for u using the Laplace transform, as in the proof of Lemma 14. 1[ we obtain 



u(t,-,y) = ~ f e xt Ri(-(3vp (K;y))d\ 

17:1 JReA=A + l 

/ e xt R^vR° x (;y))d\ (14) 

</ReA=An + l 



2^ JRcA=A +1 

exp(Aot)(^,/3^ (-,y))^ - ^- [ e xt R^(3vR° x (;y))d\, 

Zm Jr(X -e) 



where the first term on the right hand side is due to the residue at A = Ao- The first term 
can be re-written as 

exp(A t)(V>/3, PvR° Xo (- } y))ipp(x) = 

exTp(\ t)(R° Xo (pvipp))(y)ipp(x) = -exp(\ t)ipp(y)ipp(x). 

The last equality here follows from the fact that ipp is an eigenfunction with eigenvalue 
Ao, that is 

(-A - A )^ = -Pvfy. 

In order to estimate the second term on the right hand side of ffl^|) . we note that from ([7]) 
(see also (1331)) it follows that 



\R°\{x,y)\ < A(l)\V\\^\x - y\^\exp(-V2X\y - x\)\ 
if |A|, \y — x\ > I. Thus 

ll)9^(-,l/)llc^<A(e)|y|»-'|>/A|'-texp(-V2(Ao-e)|y|) 

for y <£ K', A G T(A — e) due to the fact that Rey/X > a/Ao — e for A G T(Ao — e) and 
\y — x\ > I for x G supp(t>), y £ K' . 

Hence, using the estimate on the norm of R x : C exp — > C from Lemma 13. 2\ we obtain 

\\Ri(PvR° x (;y))\\c < A(e)\V\\l-l exp(-^2(\ - e)\y\), A G T(A - e). 

Therefore, since ReA < Ao — e for A G T(Ao — e) and the factor e xt decays exponentially 
along T(Ao — e), the C-norm of the second term on the right hand side of (fb£j) does not 
exceed A(e)exp((X — e)t — \y\yj2(\ Q — e)). The term p (t,x,y) with x G K, y <£ K', 
t > 1/2, is estimated by the same expression, possibly with a different constant A(e). 
Indeed, if t > 1/2, then 



p (t,x,y) < Aex.p(-\y-x\ 2 /2t) < Aexp((A - e)t - \y - x\y / 2(\ - e)) 



since 



\y - x\ 2 /2t + (A - e)t - \y - x\y/2(\ -e) = {\y - x\/V2~t - y/(X - e)tf > 0. 

D 
We'll need additional notations in order to describe the asymptotics of p n with n > 1. 
Let al(t,y) = ipp(y) and a 2 (t,y) = exp(— et — \y\yj2(\ Q — e)). Consider all possible 
sequences a = (<Ti, ...,er n ) with a^ G {1,2}. By U™(t,yi, ...,y n ) we denote the quantity 

n"(t,z/i,...,2/n)= sup a?(t,y 1 )-...-a° n (t,y n ). 

<r^(l,...,l) 



Let Pt : C exp — > C be the operator that maps the initial function g to the solution 
p(t, •) of equation (ITTj) . Let P®g(x) = exp(A £)('?/'/3,.6')'?/'/3(a;) and P t x = P t — P^. Lemma H~T1 
states that 

||P t 1 ||<A(£)exp((A -£)i). 

The particular form of P® then implies that 

\\Pt\\ < \\Pt\\ + \\P}\\ < A'exp(\ t). (15) 

For g G C exp and n > 2, we denote 

/■DO 

!„(#) := i?J Ao 5f = / exp(-n\ s)P s gds G C. 
Jo 

Note that 

/ exp(n\ s)P t - s gds = exp(n\ t) / exp(— n\ s)P s gds 
Jo Jo 

= exp(nA £) (/„(<?) + 0(exp(— (n — 1)A £))) as t — > oo. (16) 

The functions fi, f 2 , ... are defined inductively: /i = ^g and 

n— 1 ■ 

f» = P52 k \f n L k )M vfkfn - k) ' n - 2 - (17) 

Lemma 4.3. Let K <zM. d be a compact set. For each e G (0,r) ; the function p n satisfies 
Pn(t, x, j/i, ..., j/ n ) = exp(nA t)/n(a:)^/3(|/i) • - • ^(j/n) + Qn(t, x, yi, ..., j/ n ), (18) 

where 



sixp\q n (t,x,yi,...,y n )\ < A n (e) exp(nA t)II"(t,yi, ...,y n ) (19) 



/ort> 1/2. 



Proof. For n = 1, the relation (fl8"|) coincides with the statement of Lemma [4.21 Let us 
assume that (|18[) holds for all natural numbers up to and including n — 1. A generic 
subsequence U C Y = (y ly ..., y n ) will be written as U — (zi, ..., z\u\) and its complement 
as Y \ U — (zi, ...,z n -\u\)- By the Duhamel principle applied to the equation for p n , we 
obtain 

Pn(t,;Vl,—,Vn)= Pt-s(f3v ^ P\U\(s,-,Zh-,Z\U\)Pn-\U\(s,-,Zl,-,Zn-\U\))ds 

^° UcY,U^$ 

Pt-s(/3v Y^ ex P(\ U \ X os)f\U\(-)*Pl3(zi) ■ - ■ ^fi(z\U\)) 

UcY,U^9) 



x exp((n - |Ef|)Aos)/„_[i7[(-)^8(zi) • ... • il)p(z n -\u\))ds (20) 

+2 / Pt- S (f3v ^2 ex P(\ u \^o s )f\u\(-)i>p(zi) ■ ■■■ • i/>p{z\u\)qn-\u\(s, -,z h ...,z n -\u\))ds 
^° ucy,u^q 

+ / Pt-sWv ^ <l\U\(s,-,Zl,---,Z\U\)qn-\U\(s,-,Zl,---,Zn-\U\))ds. 

^° ucy,u^q 

The second and third integrals on the right hand side of (|2"U|) contribute only to the 
remainder term. Indeed, consider the contribution to the second integral from the term 
with a given U: 

P t _ s (/3wexp(|[/|A s)/| t /|(-)V ; /3(2 : i) • ••• • $fi(z\v\)q.n-\v\(8, ;Zi, ...,z n -\u\))ds 



< AiJ>p(z 1 )...iJ>0(z\ U \) / Pt- t (v exp(|C/|A s)/|j7|(-) 

JO 

x exp((n - \U\)\ s)Tl™- lul (s,z u ...,z n _ m ))ds 
< Aip^Zi) ■ ... •if>fi(z\ U \) / exp(Ao(t-s))exp(raA s)n"" |c/| (s,Zi,...,2 n _| f/ |))ds 



JO 

< Aexp(n\ t)i!p(z 1 )...ip l3 (zi U i)'n^~ lul (t,z 1 ,...,z n _i U i) < Aexp(nA t)n£(t,2/i, ...,y n ), 

where the first inequality follows from the inductive assumption and the second one 
from ( TloT) . The third integral on the right hand side of (1201) is estimated similarly. It 
remains to consider the first integral. It is equal to 

i()p{yi) ■ ... -i/>p{y n ) / exp(n\ s)P t ^ s (f3v ^ f\u\fn-\u\)ds 

^° UcY,U=£® 

= if>p(yi) ■ ... • ^(j/n) exp(nA t) (/„(•) + 0(exp(-(n - 1)A £))) , 

where the last equality follows from ffTBT) . Thus we obtain the main term from the right 
hand side of (TT51) plus the correction 

^(yi) • ... • ipp(y n ) exp(nA *)O(exp(-(n - l)A t)) 

for which the estimate (fT9j) holds since if)p(y\) exp(— Aot) < a^(£, yi) due to (ITU]) . □ 

We will now use Lemma 14.31 to draw conclusions about the distribution of particles 
at time t. First, let us observe that for each x the sequence (J Md ip^(y)dy) fn(%), n>l, 
serves as a sequence of moments for a random variable ^' x whose distribution is defined 
uniquely. Indeed, by the Carleman theorem, it is sufficient to check that 

10 



From (TT51) it follows that there is a constant A such that 

n — 1 



i»(0)l|c<;— rlM|c«p, ™> 2 - 



Therefore, from (fT7|) it follows that for a different constant A, 

ll/»llc<^X) pr _M| IIMIcll/n-fcllc n ^ 2 > IIMIc<a 

fc=l '^ n '■ 

From here, by induction on n it follows that ||/ n ||c < A 2n_1 n!, which in turn implies ( 12~T|) 
since n! < ((n + l)/2) n . 

Let rif X (U) be the number of particles in a domain [/ C R d , assuming that at t = 
there was one particle located at x. We will write n t ' x instead of nf' 31 (K d ). Note that 

E(^((/)r = V%fc) /... [ p k (t,x,y 1 ,...,y k )dy 1 ...dy k , (22) 

where S(n, k) is the Stirling number of the second kind (the number of ways to partition 
n elements into k nonempty subsets). Formula (1221) easily follows if we write 



B.x 



([/) = J>?'*(A, 



where U = UjAj is the partition of U into small sub-domains, and then take the limit as 
maxj diam(Aj) — > 0. 

Let £P ,X be a random variable with the moments 



Let (pfi be the density on M d given by 

<Pp{x) = ippix)/ / ipp(y)dy. 
Theorem 4.4. For each x £M. d and each domain U C M. d , 

lim (exp(-X t)n?' x (U)) = ^> x I '?*{?)& 

t^oo V J J u 



in distribution. 

Proof. By the theorem of Frechet and Shohat pi)] , it is sufficient to prove the convergence 
of the moments. By (122]) the n-th moment of exp(—X t)n^ x (U) is equal to 

E(exp(-X t)n^' x (U)) n = exp(-n\ t) y^S(n, k) I ... I p k (t,x,y x , ...,y k )dyi...dy k . 



fc=i 
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First consider the contribution to the right hand side from the term with k = n. Note 
that S(n,n) = 1. We use (TT8|) for the asymptotics of p n . The contribution from the term 
q n (t, x, yi, ..., y n ) tends to zero: 

lim exp(-n\ t) / ... / q n (t,x,yx,...,y n )dyi...dy n = 0, 

t^ 00 Ju Ju 

as follows from the definition of U™(t,yi, ...,y n ). The contribution from the main term 
gives the desired expression: 



lim exp(-riA t) / .. / exp(nXot)fn(x) / ipp(yx)...ipfi(y n )dyi...dy n = E(^> x ) n ( / <pp(x)dx) 
*-^°° Ju Ju Ju 



It remains to note that the contribution from each of the terms with k < n tends to zero. 
Indeed, it is equal to 

exp(-(n - k)X t) (exp(-kX Q t)S(n, k) / ... / p k (t,x,yx, ...,y k )dyx-dy k \ . 

The expression inside the brackets tends to a finite limit as in the case k = n, while the 
exponential factor in front of the brackets tends to zero for k < n. □ 

5 The sub-critical case 

Throughout this section we assume that d > 3 and j3 G (0,/3 cr ). We will show that the 
total number of particles n t ' x tends to a random limit when t — $■ oo. Denote the integrals 
of the correlation functions by p n (t,x): 

Pn(t,x)= / ... / p n (t,x,yx,...,yn)dyx--dy n . 

Jm. d JR d 

From ©-© and (jSJ)-© it follows that these quantities satisfy the equations 

dtpx(t, x) = -Ap x (t, x) + f3v(x)'p l {t, x), 

Pi(0,x) = l 
and for n > 1: 

1 / n_1 l 

d t p n (t, x) = - Ap n (t, x) + /3v(x) I p n (t, x) + Y^ k\( n L k)\ Pi°(t, x)Pnr-k(t, x ) 

p n (0 7 x) = 0. 

As in the previous section, in order to find the asymptotics of the total number of particles, 
we will study the asymptotics of p n {t, x) as r — > oo. 
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Lemma 5.1. For (3 G (0, (3 CT ), there is a constant A such that the solution p of 177]) can 
be estimated as follows 

\\p(t,.)\\ c <A(l + t)- d / 2 \\g\\ Ccxp . (23) 

// the initial condition g G C exp in / T77]) is replaced by g = 1, then 

\imp(t,x) = l + (p8(x) (24) 

t— »oo 

in C(M d ) ; where <p fi = -R%(/3v). 

Proof. Let d > 3 be odd. As in (ITS]) , we represent the solution as an integral 



^* jReA=l 



Since the integrand is analytic in C with appropriate decay at infinity and has a limit as 
A — > 0, we can replace the contour of integration by 

7 = {z G C : Rez = — |Imz|}. 

Using the representation (JHj) for i?^, we obtain 

Pit, •) = ~ J ^{Q^ 1 + J*(A) + 0(\\\i))gd\. 

The contribution to the integral from the first term is estimated as follows 

HOJk I e^Xi-'dXWc = \\Q P d gri j ' e'a^dsWc < Ari\\g\\ Ce ^ 

where we used the change of variable s = At. The contribution from the second term is 
equal to zero since the integrand is analytic and the contour can be moved arbitrarily far 
to the left along the real axis. The third term can be treated in the same way as the first 
one, resulting in the decay in t of order t~2 _1 . The obtained estimates imply (|23|) for 
t > 1. Clearly ([23]) holds for t < 1. 

The case when d > 4 is even is treated similarly. The slight difference is that the 
contribution from the main term is now, up to a multiplicative constant, equal to 

Q p d g f e xt \-2- 1 \n(l/X)d\. 

J-y 

After the change of variable s = At, the integral is seen to be equal to 

ff I e s si~ 1 \n(t/s)ds = ri\\it I e s s^ds + t~i I e 8 s^~ 1 \n{l/s)ds. 
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The first integral on the right hand side is equal to zero since the integrand is an analytic 
function, and the contour can therefore be moved arbitrarily far to the left along the real 
axis. The second term on the right hand side has the desired order in t. 

It remains to prove (BID - Note that w(t, x) = p(t, x) — 1 is the solution of the problem 

ow(t X) 1 

— ^-^ = -Aw(t, x) + (3v(x)w{t, x) + fiv(x), w(0, x) = 0. 

By the Duhamel formula, 

w{t, ■) = — [ [ e x{t ~ s) R p x {Pv)d\ds = 

since in the domain 7 + to the right of the contour 7, the operator R x is analytic and 
decays as |A| _1 at infinity. We make the change of variables s = Xt and obtain, as t — )■ 00, 

D 

For g e C exp , we define J(g) = J Q °° P s gds. From Lemma [5TT1 it follows that 

\\J(g)\\c<A\\g\\ Ccxp (25) 

for some constant A. Let us define the sequence of functions f 1 ,f 2 ,... inductively via: 
fi = l + <pp and 

n—1 1 

f" = (3 J2 k u i kV J (vfkfn-k), n>2. (26) 

fe=i '^ n '■ 

Lemma 5.2. For eac/i n > 1 we have 



hm p n (t,-) = f n (27) 

t— >oo 



m C(R a 



Proof. For n = 1, the statement coincides with the second part of of Lemma 15.11 Let 
us assume that (j2"Tj) holds for all natural numbers up to and including n — 1. Define the 
functions Ck(t,x) = p k (t,x) — fk{x). Thus \\Ck(t, -)\\c — > as t — > 00 for k < n — 1. By 
the Duhamel principle applied to the equation for p n , we obtain 



n-l 



f n\ 

\(t, ■)= Ptsiflv Y^ U7Z ' ,.\M S > -)Pn~k(s, -))ds 



k\{n-k)V 
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ft n ~ 1 I 

= L p ^^w^w l ^ )ds+ (28) 

2 / Pt-aiPvY] jjj — '-—c k (s,-)f n - k )ds+ / p t- s (/3vy2— — '——Ck{s,-)c n -k(s,-))ds. 
Jo k=l K -\ n ~ fc J ! -/o k=l K -\ n ~ fc J ! 

From Lemma 15.11 it follows that ||P f g||c < A(l + ^) _d ^ 2 ||fl , ||c B xp f° r some constant A. 
Therefore the C-norm of the sum of the last two term on the right hand side of ( |28|) is 

estimated from above by 

ft 

-f(s)(l+t-s)- d ^ 2 ds, 

'0 

where j(s) is some function such that j(s) — > as s — > oo. The latter integral tends 
to zero when t — )■ oo since d > 3. It remains to note that the first term on the right 
hand side of (|28|) tends to f n in the C-norm, as immediately follows from the definition 
of J. □ 

Theorem 5.3. For each x G M d , the total number of particles n t ,x converges almost surely, 
as t — >■ oo ; to a random variable C with the moments m n (x) = Y^k=i ^( n ' k)fk(x)- 

Proof. Since n t ' x is monotonically increasing in t, and the moments converge to m n (x) 
(as follows from (|22|) and Lemma l5.2p . we have convergence almost surely to a random 
variable with the specified moments. □ 

6 The critical case 

Throughout this section we assume that d > 3. We will show that for (3 = /3 cr the total 
number of particles tends almost surely to a finite random limit, while the expectation of 
the number of particles tends to infinity. 

It is clear that the random variables n t ' x can be realized on a common probability 
space in such a way that n^' x < n v ' x whenever (3 < (3' and t < t! . Therefore, in order to 
show that Enf cr,x — > oo as t — > oo, it is sufficient to show that 

lim lim Enf ,IE = oo. 

This follows from Lemma 7.3 of [5], which can be stated as follows: 
Lemma 6.1. There are positive constant hd, d > 3, such that 

hmEn^ = — b J—^J x ) + 0(l) as [3 t Per 

t^ao fj cr - p 

is valid in C(M. d ), where ijjp cv is the positive ground state for L@ cr . 
Now we prove the main theorem of this section. 
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Theorem 6.2. For each x, the limit if := lim^oo n^ 1 ,x is finite almost surely. 

Proof. Let S(x) = P(r) x = oo). Let us show that S(x) depends continuously on x. Indeed, 
let x G M, d and e > be fixed. There exist 5 > and t > such that for \y — x\ < 5 the 
branching processes starting at x and y have the following properties: 

(a) The probability that at least one branching occurs on the interval [0,£] for either 
of the processes does not exceed e/3; 

(b) There is a positive function q(z) such that f Rd q(z) > 1 — e/3, which serves as a 
lower bound for both of the heat kernels p(t, x, z) and p(t, y, z). 

This shows that the branching processes starting at x and y can be coupled on an 
event of probability at least 1 — £, and therefore \S(x) — S(y)\ < e, proving the continuity. 

Suppose that S(x) is not identically equal to zero. Then S(x) ^Oona set of positive 
measure, due to the continuity. By the Markov property, S(x) > J Rd p(l,x,z)S(z)dz, 
which shows that S(x) > for all x G M. d . Let 

m = min S(x) > 0. 

xSsupp(i;) 

Let Q x = {r] x = oo}. Take iV = [4/m] + 1. Let r x be the probability that there are 
at least iV particles on the support of v before time t — 1. Clearly, there is a positive 
constant r such that 

r x > r for all x G supp(f). (29) 

Define the random times rf , rf , ... to be the consecutive instances of branching for the 
process starting at x. There is at least one particle on supp(f ) at each of the times t x , and 
t x — > oo almost surely on Q x . Therefore, from ( 1291) and the Markov property it follows 
that almost surely on Q x there is a random time t x such that there are N particles on the 
support of v at the time t x . Therefore, there exists T < oo such that P(t x < T) > m/2 
for x G supp(w). Note that T can be taken to be independent of x using the continuity 
in x of the probabilities under consideration and the compactness of supp(w ). 

Now fix an arbitrary x G supp(f). We saw that with probability at least m/2 there 
are at least N particles at time t x and therefore also at time T, that is En^"' 1 > Nm/2. 
Applying the Markov property with respect to the stopping time t x and using the fact 
that the particles move independently, we see that En^ 1 ' 1 > (Nm/2) 2 and, continuing by 
induction, that En^' x > (Nm/2) k > 2 k . 

Therefore, the expectation of the total number of particles grows at least exponentially 
with some exponent 7 > 0. On the other hand, from the arguments in Theorem 14.41 it 
follows that for /3 > /3 cr the expectation of the number of particles grows exponentially 
with the exponent Xq(/3). Since \q(/3) I as (5 ! (5 CT and nf' x depends monotonically on 
/3, we conclude that 7 = 0. Thus we come to a contradiction with our assumption that 
S(x) is not identically zero. □ 



L t 
on a common probability space in such a way that C, 13 ^ from Theorem 15.31 converge almost 



Remark 1. It is not difficult to see that the random variables u^ can be realized 
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surely, as (3 "| Per, to rf. 

Remark 2. Using the spectral techniques similar to those employed above and in [5], 
one can get the asymptotics of the higher order moments of n^' x . 

7 Appendix 

Here we prove several statements on the analytic properties of the resolvent R x . We 
largely follow [5]. The new steps concern the asymptotics of the resolvent as A — > 0. 
Denote 

A x = v(x)R° x : C exp (R d ) -► C cxp (R d ). 

The following lemma is similar to Lemmas 5.1 and 5.2 of [5] (see also p2] for a similar 
statement for general elliptic operators), the difference being that we now obtain a more 
precise asymptotics of R^ and A\ near the origin. 

Lemma 7.1. Consider the operators R\ : C exp (M d ) -^ C(M. d ) and A x : C cxp (M d ) -)• 

Cexp(J£ ). 

(1) The operators R° x and A\ are analytic in A G C . 

(2) The operator A\ is compact for A G C. 

(^) For eac/i £ > ; we /iai>e max(||i?°||, ||A A ||) = 0(1/|A|) as A -» oo ; | arg A | < ir — e. 
(4) The operator R^ has the following asymptotic behavior as A — > 0, A G C: 

R x = Q d \l~ 1 + Pd(\) + O(\\\l), d>3, d-odd, (30) 

^A = Q^i- 1 ln(l/A) + P d (A) + 0(|A|^|lnA|), d > 4, d - even, (31) 

where Pd are polynomials with coefficients which are bounded operators and 

Qdf = qd f(x)dx, 

where qd ^ 0. 

Remark. The term with the coefficient Qd is the main non-analytic term of the expan- 
sion as A — > 0. 

Proof. Let d be odd. From (JTj) it follows that the kernel R\(x) is an analytic function of 
A G C and the following estimates holds 

\R°x(x)\ < C d \\^X\ d - 2 \V\x\ 2 - d , \VJ&\ < 1, (32) 

\R°x(x)\ < C^v^V^^'llv^l 1 ^, \Vkc\ > 1. (33) 
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Moreover, these estimate admits differentiation in A and x resulting in 



,dR° 



x) 



d\ 



< C' d \V\\ d - A \VXx\ 2 - d , \V\x\ < 1, (34) 



,^( 



and 



dx ,<Ci|VAr a |e-^-l||A/Ax|^(^ + -M-), \V\x\>l, (35) 

\V x R° x (x)\ < C^\V\\ d - 2 \V\x\ 2 - d ^-, \V\x\ < 1, (36) 

\x\ 

\V x R° x {x)\ < C^VAl^le-^^llVAxI 1 ^^ + |V^|), \V\x\ > 1. (37) 

\x\ 

where C d , C' d and C' d ' are positive constants. The estimates (!32|) -( l33|) and ([34]) -( 135]) imply 
the analyticity of the operators i?° and A x . 

The estimates fl32])-([33]) and flMD-flSTD imply that the operator R° x : C exp (R d ) -+ C l (R d ) 
is bounded. Then the standard Sobolev embedding theorem implies the compactness of 
the operator A\ = vR° x in the space Cex P (K d )- 

In order to prove the third statement of the lemma, we observe that the norm of i?° 
can be estimated by J Rd \R° x (x)\dx, which is of order 0(1/|A|) as A — )• oo, | arg A | < n — e, 
due to (|32|)-(133|). 

To prove the fourth statement, we need a more detailed asymptotics of the kernel 
R° x (x) when |vAx| J, 0. Namely, it follows from the properties of the Hankel functions 
and ([7]) that there are a constant a d G C and a polynomial b d with complex coefficients 
such that 

Rl(x) = \x\ 2 - d (a d {V\\x\) d - 2 + b d {X\x\ 2 ) + 0{\VXx\ d )') as |VAx|;0, A G C. (38) 

Combined with ( )33|) and the definition of the space C cxp (M d ), this easily implies the fourth 
statement of the lemma. 

The case of even d is similar. The main difference concerns formulas (|32|) and (1381) . 
The estimate (132]) remains valid except the case d = 2, where it is replaced by 

\R° x (x)\<C 2 \\nV\x\, \Vkc\<l, 
while (138]) is replaced by 

R° x (x) = \x\ 2 ~ d (a d (V\\x\) d - 2 ln(VA>|) + b d (\\x\ 2 ) + 0(\V\x\ d \ ln{V\\x\)\)) (39) 



as |vAx| I 0, A G C. The rest of the arguments proceed as earlier, but now employing 
(159"]) instead of fl38]). □ 

The following lemma is simply a resolvent identity. 
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Lemma 7.2. For A G C , we have the following relation between the meromorphic 
operator-valued functions 

4 = R°x- R°x(I + Pv{x)Rl)-\Pv{x)Rl) (40) 

Remark. Note that (140]) can be written as 

Ri = R° x - R° X (I + PA x )-\Pv{x)Rl). 

From here it also follows that 

Rl = R° x (I + PA x )-\ (41) 

which should be understood as an identity between meromorphic in A operators acting 
from C exp (R d ) to C(R d ). 

The kernels of the operators / + /3A X , A G C, are described by the following lemma. 

Lemma 7.3. (1) The operator-valued function (I + fiA^ 1 is meromorphic in C . It has 
a pole at A G C i/ and on/y z/ A is an eigenvalue of L 13 . These poles are of the first order. 

(2) Let Aj(/3) be a positive eigenvalue of L 13 . There is a one-to-one correspondence 
between the kernel of the operator I + /3A\ i and the eigenspace of the operator L 13 cor- 
responding to the eigenvalue Aj. Namely, if (I + PA\ t )h = 0, then ip = —R x h is an 
eigenf unction of L 13 and h = flvip. 

(3) If d > 3, there is a one-to-one correspondence between the kernel of the operator 
I + (3Aq and solution space of the problem 

1 Bib 

L p (ip) = -AtP + (3v(xU = 0, if)(x) = 0(\x\ 2 - d ), -^-(x) = 0(\x\ 1 - d ) as r = \x\^oo. 

2 or 

(42) 

Namely, if (I + (3A )h = for h G C cxp {R d ), then %b = —R®h is a solution of $2]) and 

h = (3vip. 

Proof. The operator A\, A G C, is analytic, compact, and tends to zero as A — > +oo by 
Lemma [7. II Therefore (I + PAx)' 1 is meromorphic by the Analytic Fredholm Theorem. 

If A G C is a pole of (I + (3A\)~ l , then it is also a pole of the same order of R x 
as follows from ( J4T1) since the kernel of R x is trivial. Therefore the pole is simple and 
coincides with one of the eigenvalues A«. Note that A is a pole of (/ + /3Aa) _1 if and only 
if the kernel of / + (3A\ is non-trivial. Let h G C CX p(M d ) be such that | \h\ \c m imA) 7^ and 
(J + (3vR° x )h = 0. Then ip := -R° x h G L 2 (R d ) and (|A - A + [3v)ip = 0, that is ip is an 
eigenfunction of //. 

Conversely, let ib G L 2 (M. d ) be an eigenfunction corresponding to an eigenvalue A«, 
that is 

(~A-\ i )ib;+f3viP = 0. (43) 
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Denote h = (3vip. Then (|A — \i)ip = —h. Thus ip = —R\,h and ( 143|) implies that 
/i satisfies (/ + (3vR\)h = 0. Note that h G C°°(R ), h vanishes outside supp(f), and 
therefore belongs to the kernel of / + (3A\ t . This completes the proof of the first two 
statements. 

Similar arguments can be used to prove the last statement. If h e C exp (lR d ) is such 
that ||ft-||coxp(R d ) 7^ and (/ + (3A )h = 0, then h has compact support and the integral 
operator R^ can be applied to h. It is clear that ip := —R^h satisfies f )42|) . 

In order to prove that any solution of ( 142]) corresponds to an eigenvector of / + (3A , 
one only needs to show that the solution ip of the problem ( ]4"2l) can be represented in the 
form ip = —R®h with h = (3vip. The latter follows from the Green formula 

ij(x) = -(R° ({3viP))(x) + / [R° (x - yW r {y) - — R°(x - y)Hy)}ds, \x\ < a, 

J\y\=a Or 

after passing to the limit as a — > oo. D 



Remark. The relations ( H2l) are an analogue of the eigenvalue problem for zero eigen- 
value and the eigenfunction ip which does not necessarily belong to L 2 (M d ). We shall call 
a non-zero solution of fj42|) a ground state. 



Due to the monotonicity and continuity of A = Ao(/3) for (3 > (3 cr , we can define the 
inverse function 

/3 = /3(A):[0,oo)^[/3 cr ,oo). 

We'll prove that the operator — A\, A > 0, has a non-negative kernel and has a positive 
simple eigenvalue such that all the other eigenvalues are smaller in absolute value. Such 
an eigenvalue is called the principal eigenvalue. 

Lemma 7.4. The operator —A\, A > 0, has the principal eigenvalue. This eigenvalue 
is equal to 1//3(A) and the corresponding eigenfunction can be taken to be positive in the 
interior o/supp(f) and equal to zero outside o/supp(f). If d>3, then the same is true 
for the operator —Aq (in particular, j3 cr > 0). 

Proof. The maximum principle for the operator (|A — A), A > 0, implies that the kernel 
of the operator i?°, A > 0, is negative. Thus, by ([7]), for all y the kernel of — A\ is 
positive when x is in the interior of supp(t> ) and zero otherwise. Thus — A\, A > 0, has 
the principal eigenvalue (see [H]). On the other hand, by Lemma [7T3l 1//3(A) is a positive 
eigenvalue of — A\. Note that this is the largest positive eigenvalue of — A\. Indeed, if 
\i = 1/(3' > 1//3(A) is an eigenvalue of — A\, then A is one of the eigenvalues Aj of L 13 by 
Lemma [7.31 Therefore, Aj(/3') = Ao(/3) for (3' < (3. This contradicts the monotonicity of 
Ao(/9). Hence the statement of the lemma concerning the case A > holds. 

For d > 3, the kernel of — A is equal to vPd and has the same properties as the kernel 
of — A\, A > 0. Thus — A has the principal eigenvalue. Since A\ — > A as A ! 0, the 
principal eigenvalue 1//3(A) converges to the principal eigenvalue /x < oo of — A Q . On 

20 



the other hand, /3(A) is a continuous function, and therefore // = l/f3 cr , which proves the 
statement concerning the case A = 0. □ 

Remark. Let d > 3. Lemmas 17.31 and 17.41 imply that the ground state of the operator 
L 13 for (3 = (3 cr (defined by (l4"2jl ) is defined uniquely up to a multiplicative constant and 
corresponds to the principal eigenvalue of A Q . If (3 < (3 cr , then the ground state (with 
A = 0) does not exist and the operator / + (3Aq has a bounded inverse. 

We can finally proceed with the proof of Lemma 13.21 

Proof of Lemma \3. E The analytic properties of R x follow from (J4T|) and the corresponding 
properties of (I + fiAx)^ 1 which are in turn due to Lemma [7.31 

By Lemma I7.1[ the norm of A \ decays at infinity when A — > oo, |argA| < it — e. 
Therefore there is A > such that the operator (I + j3A\)~ l is bounded for |argA| < tt — e, 
|A| > A. The decay of the norm of R x now follows from ( 14 ip and the third part of 
Lemma 17.11 

Now let us prove ([8]). Let d > 3 be odd. We use ( 130]) to represent I + (3A\ as 

I + (3A X = B + C x , 

where B = I + /3A and C x = ^(QdA^ + P^A) -P d (0) + O(\\\$)). Since B is invertible 
(see the Remark following Lemma 17. 4p and C\ — > as A — > 0, A 6 C, we have 

(J + pAx)- 1 = (B + C,)- 1 = (/ - {B- l C x ) + (B-'C,) 2 - ..^B- 1 

for all sufficiently small A G C. Combining this with (J4T|) and using (1301) . we obtain (jSJ). 
In order to prove ([9]), we can repeat the same arguments, starting with ( |3T|) instead 
of (1301). D 
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